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Using tight binding band picture for 2D graphite, and the Hubbard interaction, recently we 
obtained a gapless, neutral spin-1 collective mode in graphite 3]. In this paper we present a detailed 
RPA analysis of the Neutron Scattering cross section for this collective mode. Near K— point and 
very close to F— point, the intensity of neutron scattering peaks vanishes as q 3 . This is shown using 
a simple Dirac cone model for the graphite band structure, which captures the small— q behavior 
of the system. As we move away from the F— and K— points in the Brillouin zone of the collective 
mode momenta, we can identify our collective mode quanta with spin triplet excitons with the 
spatial extent of the order of a few to a couple of lattice parameter a, with more or less anisotropic 
character, which differs from point to point. We also demonstrate that the inclusion of the long 
range tail of the Coulomb interaction in real graphite, does not affect our spin-1 collective mode 
qualitatively. This collective mode could be probed at different energy scales by thermal, hot and 
epithermal neutron scattering experiments. However, the smallness of the calculated scattering 
intensity, arising from a reduced form factor of carbon 2p z orbital makes the detection challenging. 



INTRODUCTION 

Graphite is a broad band tight-binding system com- 
posed of hexagonal graphene sheets held together by van 
der Waals interaction. This pure carbon system, inspite 
of a simplicity in electronic and crystal structure, has a 
rich physics and continues to surprise us[l]|. The four 
valence of each carbon atom is used to form a and tt 
bonds with its three neighbors. Each layer is like a giant 
molecule with resonating n bonds among many valence- 
bond configurations 0] • This structure is similar to many 
interesting systems such as carbon nano-tubes, bucky- 
ball, MgBi etc. Within each graphene layer, p z elec- 
trons of each carbon atom can hop from a site to another 
site. These electrons are responsible for the formation 
of 7r bands which touch each other at the corners of the 
Brillouin zone (Figure |2J) 

An important question is to understand the nature of 
low energy collective excitations in our zero gap planar 
p — 7r bonded graphite. We may approach this question 
in an unconventional fashion starting from organic chem- 
istry ! Planar p — it bonded molecules form the basis of 
organic chemistry, benzene being a first member. Ben- 
zene is a mini-graphite in some respects. One can also 
view graphite as an end member of planar p — -k bonded 
molecules - benzene, naphthalene, anthracene, coronene 
etc. It is well knownQ that benzene and the above se- 
quence have a spin triplet state as their first excited state. 
The next excited state is a singlet state nearly 2 eV above 
the triplet state for benzene. This remarkable singlet- 
triplet splitting, which is missing in simple Huckel the- 
ory, is a well known effect of coulomb correlation in p — tt 
bonded systems. 



A natural question is what happens to the triplet and 
singlet excitons in graphite, the end member of the above 
sequence. In a recent paper we showed that in graphite 
the triplet excitons form a well defined band in the entire 
Brillouin zone. We view the low energy part of the above 
exciton as a spin-1 collective mode in view of its gapless 
character. The singlet excitons, on the other hand form 
the well known 7r plasmon (energy ~ 7 eV) with a finite 
gap in the spectrum. Since graphite is also viewed as 
a semi-metal, our neutral spin-1 collective mode can be 
also viewed as Landau's 'Spin-1 Zero Sound'(SZS). 

Another way to understand the low energy collective 
mode in graphite is to go back to Pauling's resonating va- 
lence bond (RVB) theory of graphite^. Pauling assumed 
a dominant near neighbor singlet p — 7r bonds and their 
resonance to develop a theory of graphite. We know from 
recent developments in RVB theory that such a well de- 
veloped singlet correlation leads to a quantum spin liquid 
state. A quantum spin liquid state either has a collective 
spin-1 branch or contain spin- half spinon excitations by 
quantum number fractionization. Our finding of a gap- 
less spin-1 mode qualifies graphite to be a quantum spin 
liquid state. The gapless character of the spin-1 branch, 
however, makes it a long range RVB state rather than a 
short range RVB state as envisaged by Pauling. 

Yet another physical picture of our spin-1 collective 
mode is as follows. In graphite the valence band is 
completely filled and the conduction band is completely 
empty. The nature of particle- hole continuum in this sys- 
tem is such that a window below the particle-hole contin- 
uum opens (see figure 0}. A Hubbard type on site repul- 
sion between the electrons is in the spin singlet channel. 
This translates into an attraction in particle-hole spin 
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triplet channel. This can bind an electron and a hole 
into a spin-1 exciton; this appears as a pole for the spin 
susceptibility inside the window which exists below the 
particle-hole continuum. There are special band struc- 
ture reasons as to why our spin-1 excitations are gap- 
less, which will be explained in detail below. Use of the 




FIG. 1: (left) Lattice structure of a two-dimensional graphene 
sheets. It is a bipartite lattice composed of two sub-lattices. 
Each sub-lattice is denoted by a different color. The basis vec- 
tors of real space lattice are given by Si = -^w^eix + %e y , S2 = 



V3a; 



2 ^ x ' 2' 

§e H , where is nearest neighbor C-C distance. 
The unit cell is denoted by a dotted hexagon and contains 
two non-equivalent carbon atoms, belonging to two sub- 
lattices, (right) Corresponding reciprocal space is defined by 
b x = ^(-^ + e y ), b 2 = ^(-^| - e. y ). One possible choice for 
Brillouin zone (BZ) is a hexagonal region of the above figure. 



short range interaction for the spin phenomena is justi- 
fied, since (as in the following we will explicitly demon- 
strate for the case of spin-1 collective mode) inclusion of 
the long range tail of the interaction does not affect the 
spin physics qualitatively. 

Gapless spin-1 excitation (Goldstone modes) are char- 
acteristic of magnetically ordered systems. In the now 
popular case of La,2CuOn, a 2D quantum antiferro- 
magnet, the spin-wave bandwidth is of the order of 
~ 100 meV . But in the case of graphite, which has no 
long range magnetic order, this gapless spin-1 collective 
branch (a non-Goldstone mode) has a wide dispersion 
~ — 2 eV (figure [SJ). This energy range should be com- 
pared with triplet exciton of buckyball at ~ 1.64 eU||. 
Spin-1 collective mode of graphite exists everywhere in 
BZ and since it does not enter the particle-hole contin- 
uum, is protected form Landau damping (decaying into 
particle- hole pairs). Therefore the spin-1 collective mode 
of graphite is long-lived and could be exploited in coher- 
ent transport of spin-only (neutral) currents through the 
nano-tubes. 

Since the bandwidth of this collective mode is very 
large, each energy range should be investigated by differ- 
ent probe. Here we will concentrate on inelastic neutron 
scattering experiments. In real graphite, there are very 
small electron-hole pockets around the K— points arising 
form inter-layer hopping of the order of ~ 10 — 20 meV. 



Above this energy range, the collective mode is well iso- 
lated from the typical phonon energies and therefore the 
neutrons need not to be polarized. The low energy parts 
(regions close to T-point in figure [3J of the collective 
mode up to ~ 100 meV can be probed with thermal neu- 
trons. Hot neutrons can concentrate on higher energies 
~ 100 — 500 meV (around the T-point). The highest en- 
ergy parts ~ 0.5 eV— up- wards (regions midway between 
r — M and midway between T — K) correspond to bound 
state wave function with spatial extent of the order of a 
few unit cell (~ 2 A) and can be studied by hot and/or 
epithermal neutrons. 

However, an important practical difficulty is the small 
value of scattering cross section. As we will see in de- 
tail, the large size of the spin carrying 2p z orbital of 
carbon reduces the neutron scattering form factor signif- 
icantly, making S^q, oj) rather small in relevant regions 
in (q, uj) space. A most appropriate region, where exper- 
iment have better chance for discovering the spin-1 mode 
is in the region midway between the T and M points in 
the Brilluouin zone. 

In this paper we present the detailed RPA analysis of 
the spin susceptibility and the neutron scattering cross 
section. The bare value of Hubbard U in graphite is 
~ 8 eV, but the renormalized value for stability issues, 
should be less than 2.23i ~ 5.8 eV. We will keep U as a 
parameter to be determined by experiment and will re- 
port the results of cross section calculation for U ~ 5 eV. 
The organization of the paper is as follows. To be self- 
contained, we begin with a brief review of the cross sec- 
tion of a neutron scattering experiment which is followed 
by a summaries of single particle band picture of the 
graphene which is essential for understanding the na- 
ture of the window below the particle-hole continuum. 
The RPA formulation of the next section is applied to 
tight-binding bands of graphene in later sections. To 
get an analytical handle, we exactly obtain the spin sus- 
ceptibility within the RPA for model of a single Dime 
cone. This model captures the behavior near the T and 
A'— point. We calculate the intensities of the neutron 
scattering peaks for the linearized model and also discuss 
the effect of including long range tail of the interaction. 
Finally we report the numerical calculation of the inten- 
sities of the neutron scattering peaks for entire BZ along 
with calculation of the real space profile of the excitonic 
wave-function. To be self-contained, in appendix A, we 
review excitonic formulation of the problem. 



NEUTRON SCATTERING CROSS SECTION 
BASIC FORMULATION 

Neutron scattering is a very weak probe and does not 
disturb the target very much. It makes the linear re- 
sponse theory a very useful tool to interpret the neu- 
tron scattering data. Therefore according to fluctuation- 
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dissipation theorem, the interpretation of neutron scat- 
tering experiment requires, only an understanding of the 
target itself. In this section following ^5| , we briefly sum- 
marize the formulation of interaction of neutrons with 
matter. 

The wave-length in 1 and energy in meV of neutrons 
are related by 



A = 9.04 £T 1/2 



(1) 



Low energy neutrons are often described as cold, ther- 
mal, hot, or epithcrmal. The typical energy scales and 
corresponding de Broglie wave-lengths are given below 





E/meV 


A/1 


Cold 


1 - 10 


9-3 


Thermal 


10 - 100 


3 - 0.9 


Hot 


100 - 500 


0.9-0.4 


Epithermal 


500- 


0.4- 



(2) 



A simple golden rule for the scattering of neutrons 
gives the cross section proportional to dynamical form 
factor 5(q, w), which is itself proportional to Inrx(q, oS). 
The dynamical susceptibility Imx(q, u), measured by 
neutron scattering yields incisive information about the 
electron-electron interactions on precisely the relevant 
energy and length scales. It is widely considered as one 
of the most important probes of strongly correlated ma- 
terials. 

Assume that the incident neutron gets scattered from 
initial plane- wave state |k) to |k'), whereby the target 
acquires a transition form initial state \n) to final state 
|n'). Within the Born approximation, the contribution 
form this process to total cross-section is given by |l5j 



/ d 2 a \ _ V 

[dndE'J^ k 



k'n'\V\kn) 5(hu - E n , n ) (3) 



where E n i n — E n i — E n is the energy change of target, 
and 



(k'\V\k) = ^ [ dr e *( k - k ')- r V(r) 
2ttH J 



(4) 



For none zero temperatures, there is a range of accessi- 
ble initial states with probability p n ~ e~ l3En . The ba- 
sic quantity which is measured, is the partial differential 
cross section 



d 2 a 
dQ.dE' 



jY,P»\ ( k '«'I^IM| 2 K*»> ~ E n'n) (5) 



In the above expression V is the interaction potential 
of neutrons with the target. The neutrons may inter- 
act with ions either through nuclear forces of spin and 
orbital angular momentum dipole interactions, or, with 
electrons dipole. The case of relevant to our problem will 
be interaction with the spin dipole moment of 2p 'band 
electrons'. 



Scattering by electrons 

If we forget about the orbital motion of electrons, the 
magnetic field of its spin dipole moment /i e = — 2/iss 
becomes 



H = V x I /x e x 



R 



(6) 



so that the interaction potential between the electron and 
neutron becomes 

- XT ~ v-7 s x R, 

- IfJ-N o-.H = 2j^ N ^ B c. V x (7) 

where a denotes the Pauli spin operators. Plugging it 
into equation JSJ and after some algebra 01 we nn d 



dfldE' 



m N 
2ttH 2 



(2-/fJ, N fi B ) 2 ('in) 2 — 22 PnPa 



{na\(a.QY\n'(r')(n'(T'\a.Q\na) 8{hu - E n , n ) (8) 
where the operator Q is defined by 

Q = ^e«qx(s,xq) (9) 

i 

Here q = k' — k is the momentum transfer, and q = q/|q| 
is the unit vector along q. Finally putting every thing 
together, we have 



dttdE' 



k' 

rg-5(q,w) 



(10) 



where ro = je 2 /m e c 2 and the structure factor <S(q, u>) is 
given by 

a/3 nn' 

The coefficients f a p depend on whether we are using po- 
larized or unpolarized neutron beam and are given by 



fa(3 = ^2p a (cr\a a ap\<j) 



(12) 



For experiments with unpolarized neutrons f a p = 5 a p 
while for polarized neutrons, say with polarization | f), 

( I i \ / 1 \ / i \ 
/«?= -U0 = 010 I -tOO (13) 
\0 01/ \ 1 J \ / 

where the second term really distinguishes polarized ex- 
periments form unpolarized (Sap) one. In polarized neu- 
tron scattering investigation of spin-1 excitations, one 
counts only the neutrons that after scattering, have 
flipped their spin. However, fortunately in case of 
graphite the spin collective mode exists up to electronic 
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energies of the order of eV, where there are no phonon 
degrees of freedom (phonons inhabit in energy scale tens 
of meV) with comparable energies. So one can safely do 
the experiments with unpolarized neutrons. With these 
assumption e quation can be written in terms of the 
Fourier component of the spin density operator as 



EOW-fe^)^ / dt e-^iS^SUt)) (14) 

a/3 



where in the last step we have used the Fourier represen- 
tation of the Dirac delta function in equation and 
the interaction picture representation of the spin opera- 
tor, that is: 



-?q.r 



(15) 



Itinerant electrons 

For the band electrons, S^ a is represented in Bloch 
basis as 

J kk'AA'ss' 

(16) 

where A is the band index, s =t, J, is the spin projection 
and ip£ represents the Bloch wave-function. Using equa- 
tions 1|32|1 and after some algebraic steps, the repre- 
sentation of the spin operator in Bloch basis of graphene 
becomes 



S la — o y) /l 



AA' Af -a A' 
k,q ^k+qs u ss'^Ws' 



(17) 



kAA'ss' 



where (A, A' = ±1) 



/k,q - 2 F ° (Cl 

A' 



1 



1 + XX' e iv>( - k+ ^~ iv>( - k ^ 



+ ^. J F 1 (q)e- l¥ ' (k)+lk - d 
-l-^ j p*(_q)g ,; ^( k +q)- i ( k +q)- d 

The form factors and phase factors are given by 

»<p(k) = /( k ) 
l/(k)| 

Fo(q) = J dr\^r)\ 2 e^ r 
Fi(q) = J dr <f (r)e iq - r 0(r-d) 



(18) 

(19) 
(20) 
(21) 



where /(k) is defined in equation (|34J) . </>(r) and </>(r — d) 
represent the atomic p z orbital at some carbon atom and 
its neighboring site, respectively. 



Neglecting the form factor F\ (q) which mixes the two 
neighboring orbitals, and using the symmetry between 
two sub-lattices of graphite, equation (|T7|) can be written 
as 



S qQ =F (q) 



kAss' 



Af l„c 
C k+qs 2 CTs 



/C ks 



(22) 



This is in fact the one-band result which is not surprising, 
because we have neglected the cross term between bands, 
and therefore the final result is same as one-band case, 
with of course summation over two separated bands. 

Using the fluctuation dissipation theorem 0, 0] , we 
have 
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(1 - e-^)5(q,u) - — |F (q)| 2 Im X (q,w 

Hence, equation (|10|) can be written in the form of 
tPa rl k' , 



(23) 



dVLdE' 



tt(1 



) k 



|fb(q)rimx(q,w) (24) 



Therefor the peak structure of neutron scattering cross 
section is reflected in dissipative part of the susceptibil- 
ity, with an overall atomic form factor, reflecting the 
band nature of electrons involved in the process. The 
Imx(q, to) includes the effects of interaction among elec- 
trons. Treating the effect of interaction in RPA ap- 
proximation, since the life-time effects are beyond RPA, 
Imx RPA has very isolated sharp peaks at resonance fre- 
quencies w s (q), that is 

Im X flPA (q, w) « Z(q) 6(hu - foj.fq)) (25) 

which defines a dimension-less quantity Z(q). By 
Kramers-Kronig relation, the above equation is equiv- 
alent to 



Re X RPA (^iv) 



1 



Z(q) 



77 H I 



w a (q)' 



w«w s (q) (26) 



Within the RPA, one does not need to perform any fitting 
to obtain Z(q). Upon using equation l|42|) . Z(q) can be 
obtained in terms of bare susceptibility as 



*(q) = ^ 



Jj ReX ° (q ' w) 



(27) 



J=w a (q) 



The measured intensity at an angle corresponding to mo- 
mentum transfer q, apart from overal factors is given by 



J(q) = |F (q)| 2 Z(q) 



(28) 



Now let us discuss the effect of atomic form factors, which 
directly affects the intensity of the neutron peaks. For 
2p z atomic orbital, we have 



z e 



-Zr/2a 



(29) 
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where do is the Bohr radius and Z is the effective nu- 
clear charge, which for 2p orbitals in carbon is w 1.56. 
Assuming that, there is no momentum transfer along the 
direction perpendicular to graphene planes, the atomic 
form factor becomes 

*b(q) = ^5-2— ~ (30) 

Note that ao ~ a/ 5 where a is the lattice parameter 
as depicted in figure ^ F° r momentum transfers much 
away from the center of the BZ, the weight of the neutron 
scattering peak is going to be reduced by |Fo(q)| 2 . 

For high energy part of the collective mode branch, we 
will need high energy incident neutron beam. Consid- 
ering the neutron kinematics, the lowest g-values which 
can be attained (at the lowest possible scattering angles 
of ~ 2 — 3 degree) are relatively high, i.e. much higher 
than the q— vector of the first Brillouin zone. Hence, 
one needs to know the magnitude of the form factor at 
q— values from ~ 5 — \2°A~ X which corresponds to 
qa ~ 12 — 29.5. The above estimate of the atomic form 
factor indicates that the intensity of peaks is essentially 
zero beyond the first BZ!. This makes the detection of the 
high energy parts a technically demanding task and the 
appropriate region of the BZ should be chosen carefully. 
Keeping this consideration in mind, after calculating the 
coefficient of the delta peaks of Imx(q, u) in RPA ap- 
proximation, we will suggest the best regions of the BZ 
(figure EJ to be probed by neutrons. 
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FIG. 2: 7T— bands of graphene along TKMT loop, at tight- 
binding approximation. The M— point is a saddle point and 
gives rise to a van-Hov singularity in single-particle DOS. In 
this figure energy is in units of t and the horizontal axis de- 
notes qa, with a defined in figure 



(33) 

with 

/(k) = e l ^T + 2e - i i5s cos(^) (34) 



THE TIGHT-BINDING BAND STRUCTURE 

Within the tight-binding approximation, and neglect- 
ing the overlap of 2p z atomic orbitals of neighboring car- 
bon atoms, the band dispersion of graphene are given by 



Hi 



c/v 



±t £k where 



£k = y 1 + 4 cos(^-) cos(-|-) + 4cos2(^-) (31) 

where + and — signs correspond to conduction and 
valence bands, respectively. Here a is the length of 
translation vector in one of the sub-lattices of graphene: 
a = \/3x (C-C bond length) as in figure The above 
bands touch each other at the if —points (corners) of the 
BZ as depicted in figure [2J Figure [3] shows the DOS of 
the dispersion (|31|l . 

For later reference, the normalized Bloch orbitals of 
graphene at tight-binding approximation become || 



= E ^XrW 



(32) 



Re A 



(p - R) ± 7^'" d *(p - R - d) 



where R runs over one of the sub-lattices of graphene and 
R — d refers to carbon atom in neighboring site which 
lives in the other sub-lattice. Here <j) is normalized p z 
atomic orbital. 

The above bands when linearized around if— points of 
the BZ become 



where Kvf 




conduction band 
valence band 



(35) 



_ V3 



ta. The DOS corresponding to the this 



linear dispersion is given by 



2ttK 2 



So the p(e) oc |e| pseudo-gap arises form the linear fea- 
tures of the band structure near the K~ points. It can 
be shown that, within DMFT, this cone-like feature re- 
mains robust against the increase in the on-site repulsion 
U |9j . Hence the on-site interaction does not destroy the 
Dirac cone picture. It only renormalizes the Fermi veloc- 
ity vf ■ Figures [21 an d El depict the it bands of graphene 
and corresponding DOS. 
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FIG. 3: The DOS corresponding to tight-binding band equa- 
tionl KUI . The energies are in units of t. The linear behav- 
ior for small energies results form the cone-like features of 
the band structure near K— points and the saddle point at 
M— point is responsible for the cusp in DOS at \uj\ = 1 



THE RPA SUSCEPTIBILITY 

In order to calculate the contributions to total RPA 
susceptibility form different channels, we need the (re- 
tarded) bare susceptibility which is given by 



,o 



(q ->= ^-(t;4 + .o+ <36) 



where N is the number of unit cells. At T = the 
conduction band is empty and we have /k+q = 0. At 
half filling which is the case for undoped graphene the 
valence band is completely filled and we have /k = 1. 
Therefore the bare susceptibility becomes (energies are 
in units of t, unless otherwise specified) 



k 

= -Y, 



-1 



N ^ hu - (e£ +q - el) + z0+ 



1 



(£k+ q + £k) -fuv- i0+ 



(37) 



The imaginary part of x° upon using the formula 

1 1 

— — = V — V inofx) 

x — iQ + x 

can be written as 

W(q, W ) = ^^J[^-(4 +q -4)] 

k 

= [ d 2 kS [Ku - (£ k+q + £k)] (38) 
47r Jbz 

where A = v ^ is the area of hexagonal unit cell. Us- 
ing Kramers-Kronig relation, the real part of %° can be 



FIG. 4: Particle-hole continuum of graphite. Horizontal axis 
is center of mass momentum of particle-hole pair and vertical 
axis is the particle-hole energy. Gray region is the particle- 
hole continuum, for the spectrum given by equation (f3ipi . The 
energies are in units of t. This is obtained for a fixed q, by 
finding the maxima and minima of the particle-hole energy 
£q+k + £k for random walker k walking in Brillouin zone. 



obtained from the imaginary part 

7T J_ oc W' - UJ 



(39) 



The particle-hole continuum is a region in lu — q space 
in which integrand of equation 1)38(1 is none-zero. Out of 
this region the denominator of summand which gives x° 
is none-zero and the integral becomes well-behaved and 
easy to be performed numerically. 

The contribution form the triplet particle-hole channel 
to the RPA susceptibility (particle-hole ladder summa- 
tion) is given by the equation 



x(q,w) 



x°(q,w) 



1 - u(q)x°(q,u;) 



(40) 



where v(q) is the Fourier component of potential. The 
short range part of the interaction is sufficient to account 
for the spin phenomena. Hence we use the standard Hub- 
bard model for this tight-binding system, 

H = -t ( c L c ^+4^) + u 12 n ^ ni i ( 41 ) 

<z,j>,cr i 

with repulsive interaction between the particles U > 0. 
For this interaction v(q) = U, and the contribution form 
triplet particle-hole channel to the RPA susceptibility be- 
comes 



x(q,w) 



x°(q,^) 

l-tfX°(q,w) 



(42) 
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The collective mode can exist if the above RPA suscep- 
tibility diverges. The roots of the denominator are the 
solutions to equations 

Im X ° (q, u) = 0, Rex (q, w) = - (43) 

One can see from equation 1)3 7|l that out of particle-hole 
continuum where hoj ^ £k+ q + £k, we have Imx° = 0, i.e. 
X° is purely real. Moreover below the particle-hole con- 
tinuum ha; < £k+q + £k and thus from equation 1)37(1 one 
can see that Rex = x° > 0. This shows that below the 
particle-hole continuum, there could be a solution to the 
collective mode equations in triplet particle-hole channel. 

Hence opening of a window below the particle- hole con- 
tinuum provides the unique opportunity for the existence 
of a solution to equation l|43l) . The only remaining part 
is to see if this solution is out of the particle-hole contin- 
uum or inside it. To this end we need to solve equation 
(|43|l numerically. Our strategy is to fix q and then for a 
given U look for a value of omega which fulfills the col- 
lective mode criterion, equation (|43|) . Figure [S] shows the 
numerical solution of the collective mode equation. 

Before going into the numerics, it is instructive to work 
out a simple model which is relevant to the low-energy 
physics of the problem. The asymptotic behavior of the 
collective mode dispersion near T and if —points can be 
understood in terms of a very simple model. The clue 
is that near the if— points of the BZ, the spectrum is 
(Dirac) cone-like and this linearity holds true up to ~ 
O.bt. So as a model it seems plausible to start with the 
simplified model of Dirac fermions, living in a circular 
BZ and subject to Hubbard interaction. It turns out 
that this problem will teach us lots about the nature of 
spin collective mode. 



see the qualitative difference between the case of 
a isolated graphene sheet and a real graphite in 
which the long-range (small— q) tail of the interac- 
tion becomes essential. We will include these effect 
analytically and will see that at T— point which is 
susceptible for instability, the collective mode still 
survives. 

Moreover, it will teach us about a mechanism by which 
the system exhibits ID characters. It turns out that a 
very same mechanism is responsible for the existence of 
the spin collective mode in almost large fraction of BZ 
area. 

Starting with equation (|38|l 

W(q,w) = [ «[n«-n»F(|k + q| + |k|)] 

47r J\k\<k, 



A 



ds 



inhvFj |V fc / q (k)| /q(k) 



(44) 



where 



/ q (k) = |k-q|+|k|, 



Vp 



and ds is the length element on / q (k) = z. The equation 
/ q (k) = |k— q| + |k| = z defines and ellipse with principal 
axis equal to | and half conical distance equal to | = 
Therefor Imx° is none zero if and only if 



z > q uj > VFq 



particle-hole continuum (45) 



This result should be compared with figure 0] around 
if —point and T— point. 

Assuming that cf> is the angle between q and k and 
k = |k|, the equation of ellipse becomes 



SIMPLIFIED MODEL OF CONE-LIKE BANDS IN 
CIRCULAR BZ 

In this section we use the linearized spectrum of equa- 
tion as a model to evaluate x° analytically. We as- 
sume that for this model BZ is a circle of radius k c which 
for the low-energy regime can be thought of as infinity. 
We assume that there is a single cone [Tjj at the center of 
the circular BZ. Note that in real graphite there are two 
such a cones. Such a simplified model seems to capture 
the essential low-energy physics of band structure (|31|l . 
The linear spectrum of equation (|35|l holds true over a 
very large range of energies (up to £ ~ 0.5t) which can 
be seen qualitatively from the density of states (figure|3J). 
Such a model will be relevant to 

• Wave vectors near if —point, of course with modi- 
fied parameters such as Fermi velocity vf , etc. 

• Wave vectors near T— point. The analytical knowl- 
edge gained about the small— q behavior of the dis- 
persion of spin collective mode will enable us to 



k = 



z 2 -q 2 



2{z — qcos4>) 



(46) 



which gives 



dk 



q(z 2 — q 2 ) sin</> 
2(z — qcos (j)) 2 



ds = \Jk 2 d§ 2 + dk 2 



z" - <r 



|V fc / q (k)| 



2(z — qcos 4>) 2 
2(z — qcos 



2zqcos( 



\J z 2 + q 2 — 2zq cos 
Putting every thing together we obtain 

z 2 — q 2 f z 2 + q 2 — 2zq cos < 



(47) 



Imx°(q,w) = 



IdirhuF 



(z — qcoscj)) 3 



2z 2 - q 2 



16hv F y/z 2 - q 2 
V3a 2 uj 2 



i6hv 2 F 



2Vf<1 , qv F < uj < (2k c + q)v F (48) 



uj 2 — vpq 2 
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where in last step we have used A = ^a 2 to restore the 
lattice parameter a. The real part is obtained using the 
Kramers and Kronig relation, eauation (|39[l . along with 
change of variables uj' = qv F coth rj, coth 770 = 1 + 
To leading order in 2 ^ £ - we obtain 

& q 



Re X °(q,^) 

V3a 2 q 2 
32-Kh 

Via 2 

2 2 

q z v F - 



V3a 2 
32fmvp 



2k " VF duj' 2uj 2 



2 2 

v F <r 



u y/uj 2 — v 2 F q 2 



2 coth 2 77-I 



sinh r\(yj — qvF coth 77) 



Ah 



2k c 



[wlog( ) - qv F { V 1) 



2uj 2 



: (arctan 



qv F 



\/(fvp 
V3k c a 2 
8nhvp 16irhvp v / q 2 v 2 7 



2uj 2 



qv F + 

- q 2 v 2 F 



-)] 

2 n 



arctan 



uj' 



qvF 
qv F 



qvp > uj 



(49) 



Now going back to equations l|43|) we see that in region 
qvp > uj where the above formula for Rex°(q, oj) is valid, 
the imaginary part is identically zero. So the dispersion 
of collective mode is the solution to equation 



1 

U 



V3k c a 2 
8irhvp 



Via 2 2uj 2 



16nhvp 



' Vp = arctan J 1 ' 



qvp — uj 
(50) 

Near the particle-hole boundary the RHS of the above 
equation can be expanded to obtain (note that Tivp — 
taVl/2) 



Vla 2 q 3 / 2 



0{q Vl 



UJ 



,1/2 



(51) 



The above equation has a solution, passing through the 
origin for any value of U < U\ — . To get some idea 
about the value of U\, let us choose k c in such a way 
that the area of our circular BZ is equal to the area of 
hexagonal BZ of the original problem, that is 



irkt = 



8tt 2 
VI c 



k c a = 2 



,1/2 



2tt 

71 



1/2 



3.81 



which gives U± = (27r\/3) ' t ~ 3.30<. Hence to leading 
order in qa the dispersion of the spin-1 collective mode 
near the T— point becomes 



(52) 



w s (q) = v F q(l - o?a 2 q 2 ) = v F q - co B (q), 
as u> — > v F q — * with 

1 UU X 

Note that equation l|52|) makes sense only if |a 2 a 2 g 2 1 <C 1. 
In particular the smallness of a implies that U must 



be below U\ and far enough from it. In fact, numer- 
ical calculations for the full band structure shows that 
for U > U c w 2.23 t, an instability emerges around the 
T— point (see figure^. Since the linear model is relevant 
to T— point (and also the if— point, but with renormal- 
ized v F ), the requirement of U < U\, warns us of an 
instability. However, since the real graphite is stable, the 
normalized value of interaction must be less than 2.23 t 
and we are not worried about it. 

Here Uujb is the binding energy of particle-hole pairs 
with center of mass momentum q. The repulsive interac- 
tion U among particles becomes attractive for particle- 
hole pairs in triplet spin state and binds them together to 
form a bound state with binding energy Hujb- To see why 
there exists a bound state for arbitrarily small attraction 
U, note that, right above the particle-hole continuum 
and very close to the continuum boundary, Im^ can be 
written as 



Imx°(q,cj) = 



V3a 2 



-3/2 



327U 



(53) 



with a square root divergence at the lower edge of the 
particle-hole continuum in u> — q space. This expres- 
sion has the same form as the one dimensional density 
of states (with energy measured from fiqv F ). Note that, 
in fact, Imx (q, w) — irpq(uj), where Pq(uj) is the free 
particlc-holc DOS for a fixed center of mass momentum 
q. That is , the particle-hole pairs have a phase space for 
scattering which is effectively one dimensional. Thus we 
have a bound state of particle-hole pairs in spin-triplet 
channel for arbitrarily small interaction U. However, we 
also have a pre factor q 3 ^ 2 which scales the density of 
states. This together with the square root divergence 
of the density of states at the bottom of the particle- 
hole continuum gives us a bound state for every q as 
q —> Q, with the binding energy vanishing as ~ q 3 as 
shown above. Equation l|53|l shows that entire mecha- 
nisms responsible for the formation of spin-1 zero sound 
(SZS) near the T and if— points, is due to a region near 
the bottom of the particle-hole continuum. In fact for 
center of mass momenta close to T, if— points, there will 
be a one dimensional manifold on which the particle-hole 
energy has its minimum and this one- dimensional mani- 
fold of minima is responsible for square root divergence in 
particlc-holc DOS. The M— point has a similar property. 
These explains why for center of mass momentum corre- 
sponding to M— point too, the spin-1 collective mode is 
there for arbitrarily small interaction U . 



neutron resonance peaks for the cone model 

Plugging equations (|49|l and l|52|l into the above equa- 
tion, it is straightforward to show that near the T and 
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K— point the peak intensity behaves like 



3^ 4 Ut 2 ( q x 3 



4 ([7i - C/) 3 V fc c 



1 - 15tt^ 



[/ 2 



. 2 / \ 3" 



(Ui -U) 2 \k c ) ' ' \ k 



(54) 



This equation shows that for a graphene layer the res- 
onance peak residue, Z(q) vanishes as ~ q 3 , where q 
is measured from T or K— points of the BZ. The atomic 
form factor is essentially constant within the first BZ and 
therefore the cross section for momentum transfers near 
r and K are very small. Vanishingly small Z(q) accord- 
ing to appendix A, means spatially large wave-function 
which can not be excited by neutrons with a typical de- 
Broglie wave-length of ~ 1A. 

Due to the presence of the other graphene layer and 
screening arising from the interlayer hopping among the 
layers, the power law Z(q) ~ q 3 may change in real 3D 
graphite. In next subsection, we will discuss the effect 
of semi-metallic screening, which amounts to taking into 
account the long range tail of the interaction among the 
electrons. 



effect of long-ranged part of the interaction 

Having established the existence of a gapless spin-1 
collective mode branch within Hubbard model and the 
RPA approximation, we will discuss whether the semi- 
metallic screened interaction of 3D stacked layers will 
affect our result. In tight binding situation like ours, 
the spin physics is mostly captured by the short range 
part of the repulsion among the electrons. So we do not 
expect a drastic change in the behavior of spin collective 
mode, if we correct the interaction, by including long 
range tails, i.e. modifications near q = 0. Below, we will 
show this explicitly and will discuss similar situations 
using an asymptotic analysis. 

Interaction, including interlayer scattering between 
layers separated by distance d is given bv(lfj| 



v(u>,q) 



2ire 2 



sinh(<7<i) 



£o<7 



[cosh(gd) + S mh(qd) X o(u>, q)] 2 - 1 



Starting from equation l|3()l) , one can see that the value of 
Rex ([j = 0,q = 0) is finite. Hence we have the leading 
behavior v(0,q) ~ q^ 1 / 2 . This behavior is something 
intermediate between the metallic screening, i.e., v(q) ~ 
(q 2 +k 2 cr )~ 1 and insulator behavior, i.e., v(q) ~ q~ l . Our 
argument below will assume the power low v(q) ~ q~ v 
for long range behavior of screened interaction. 

The essential feature of the cone model which lead to 
a gapless bound state is the square root divergence in 
DOS of free particle-hole pairs. So let us encapsulate the 



essential features into the following form 



Im X °(q, W ') 



- q 



(55) 



which accommodates key features of the problem, 
namely, square root divergence and overall q v depen- 
dence. This enables us to study the interplay between 
these two features and their effect in modification of 
graphene solution. Assuming the collective mode dis- 
persion to be 



w s (q) = q - cq 



(56) 



Using the Kramers-Kronig relation, the collective mode 
equation becomes 

1 n , , , > f dio' q 

= Re X °(q,w s (q))' ' 



v(q) 



u 2 + cql 3 



q" 



lu' - LU s (q) y/up - q 

du 



d-\fuP — q 

w'-cj s (q) q J u 2 + q-oj s (q) 
du 

q 



which gives 



(57) 



This equation is important in that, it determines the 
qualitative difference between the dispersion of the collec- 
tive mode in graphene and graphite. In our case rj = 3/2. 
The graphene is characterized by v = 0, or equivalently 
(3 = 3 which is nothing, but equation Q52[l. Once we add 
the other layers of graphene, i.e., v = 1/2, we obtain 
f3 = 2. This means that, at a given wave-vector, the 
effect of semi-metallic screening is to increase the bind- 
ing energy. Another interesting observation is that an 
insulator type of screening (y = 1) manages a collec- 
tive mode with strictly linear dispersion but with slope 
1 — c. So if c is positive, the collective mode is below 
the particle-hole continuum and survives. However, if 
c is negative, it enters the particle-hole continuum and 
decays into particle-hole pairs. 

The case of strictly one- dimensional problem corre- 
sponds to r\ = 0. So a none zero r\ is a peculiar fea- 
ture of our particular problem which should be taken 
into account in putting ID-like features on formal basis 
in order to make it appropriate for bosonization. It is 
seen from equation l|57|l that for the case of a real ID 
problem (3 = —1v. Since in semi-metallic and insulator 
case v > 0, the exponent f3 becomes negative and hence 
the expansion (|56|l makes no sense. Therefor a bare ID 
DOS (square root divergence) is not sufficient to bound 
electron- hole pairs together to form spin-1 zero sound 
(SZS). To conclude, both square root and g 3 / 2 pre-f actor 
are essential to manage a gapless spin-1 collective mode 
in graphite. 
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Now let us use a similar asymptotic analysis to see how 
the leading behavior of Z(q) gets modified by going from 
graphene to graphite. Differentiating Kramers-Kronig re- 
lation, equation (|39|) with respect to to and using (|55|) 
gives 



U)—uj s (q) 

du 



(u 2 + q-uj s (q)y 



= <?" 



dy/uj 1 — q 
(c'- Ws (q)) 2 
du 

(u 2 + cqf ) f 



Using equations l(T7|) and (|57|l we obtain 
Z(q) ~ g 2f '- 3 ^ 



7 ^-3/3/2 



(58) 



The case of graphene corresponds to n = 3/2 and v — 
which gives rise to q 3 dependence in agreement with 
equation (|54(l . This equation says that in real graphite 
(u = 1/2), Z(q) - q 3 ' 2 , i.e. the effect of semi-metallic 
screening is to strengthen the intensity of resonance peak 
at a given wave- vector. 

The above analysis reveals that, following properties 
are necessary to conspire to manage a spin-1 collective 
mode for the cone model of graphite: 

1. Dirac cone spectrum, i.e. a pseudo-gap of the form 
g(e) = |e| to allow for a window below the particle- 
hole continuum and to produce a semi-metallic 
screening v(q) ~ q~ x / 2 for 3D stack of graphite. 

2. Square root divergence of the DOS of free particle- 
holes at the edge of particle-hole continuum, ac- 
companied by q 3 / 2 measure which comes form two 
dimensionality of the original problem. 



NUMERICAL CALCULATION OF THE WEIGHT 
OF NEUTRON SCATTERING PEAKS 

So far our cone model which was relevant to small-g 
behavior of the collective mode was especially suitable to 
address the physics near the T— point (and also K— point, 
with renormalized parameters), e.g. the question of mod- 
ifications arising form inclusion of the long-range tail of 
interaction which arises form the stacks of 3D graphite. 
But for the neutron scattering experiments, we need to 
repeat the calculations in entire BZ. To find the reso- 
nance frequencies w s (q) in entire BZ, we have to numer- 
ically solve for equation l|43|) . In order to do so we have 
to perform the summation (|37|l or equivalently 



Re X °(q,w) 



V3 a 2 
8tt 2 



d 2 k 



BZ (£k+q + £k) - TlW 



(59) 




where A = ^a 2 is the area of units cell. We know form 
our analytical solutions that, below a certain value of 



no 



FIG. 5: The spin-1 collective mode for different values of U . 
Note the emergence of instability at T— point at U ~ 2.2. 
Also note the asymptotic behavior of the dispersion of the 
spin-1 collective mode near F and A"— points. According to 
analysis of the single cone model, the collective mode dis- 
persion near these points asymptotically approaches the con- 
tinuum boundary. At M— point the resonance energy is at 
~ 0.83 X t « 2.0 eV and is well isolated from the boundary 
of p-h continuum. Lower energy regions are towards T and 
if -points. Towards F-point binding energies are larger and 
better for neutron scattering as discussed in the text. 



U, the solution to collective mode equation (|43|l lies be- 
low the particle-hole continuum and asymptotically ap- 
proaches the lower boundary of continuum at if and 
T— point, Hence the denominator of the integrand in the 
above equation in region in which we must look for the so- 
lutions w s (q) is none zero and with full band dispersion 
of equation l|31|l . the integral can be done numerically, 
without any problem. However, numerically it becomes 
more and more difficult to get the square root divergences 
as we go closer and closer to if —point. This is because 
of the pre- factor q 3 ! 2 in equation (|53() . 

Figure [5] shows the numerical evaluation of the disper- 
sion of spin-1 collective mode for different values of U. 
As can be seen in figure El near if— point the collective 
mode pole is very close to the continuum boundary. One 
can fit a dispersion of type l|52l) near T and if— points, 
but the coefficients will be different from that given by 
equation l|52|l . Since near the if— point the binding en- 
ergies are very small, the particle-hole pairs are bound 
loosely and the spatial extent of their wave-function is 
large which makes it difficult to excite them by neutrons 
with typical wave-length A ~ A. Therefore we suggest the 
neutron scattering experiments, NOT to focus around the 
K— points. For the same reason the region very close to 
T— point should be avoided. 

Note that in figure^ at M— point, there is a collective 
mode for any value of U. In fact at this point, the DOS 
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FIG. 6: Plot of /(q) = \Fo(q)\ 2 Z(q) as a function of momen- 
tum transfer along YM direction. The TV distance is fa 3A~ . 
Because of the effective nuclear charge Z — 1.56 for 2p z or- 
bitals, the form factor decays very quickly, so that beyond 
the first BZ it is almost zero. This figure suggests that the 
M point is not suitable for neutron scattering. Since momen- 
tum transfers in 2 nd - 4 th BZ (5 - 12A" 1 ) are involved at 
which peak intensity is almost zero (see the discussion follow- 
ing equation IM0|l . 



FIG. 7: The approximate bound state wave-function within 
the RPA approximation. The center of mass momentum q is 
at 0.2 YM « 0.3A -1 . Unit of length is the lattice parameter 
a ~ 2AQA. Contours corresponding to |^ q (r)| > 0.1 have 
been plotted in the base. Note a very soft anisotropic pattern 
of contours. The normalization of the wave-function is such 
that V'q(O) = 1. This figure explicitly demonstrates why at 
low energies we have a collective mode. We expect ~ 0.5 eV 
neutrons to be able to excite such a bound state. 



of free particle- hole pairs diverges as ~ l/y/u) — 2t, where 
2* corresponds to bottom of the particle-hole continuum 
at M. Bindings at M— point and its neighborhood, es- 
pecially towards T, and also the point mid-way between 
TK and region towards T are strong enough to give rise 
to a particle-hole pair, with small enough wave functions 
that can be excited by thermal hot and epithermal neu- 
trons. At M— point, the resonance energy is ~ 2 eV, 
but since at these energies momentum transfers beyond 
the first BZ are involved, the effect of atomic form factor 
reduces the intensity of the delta peak almost to zero. 
This makes neutron scattering at M-point a very diffi- 
cult task. Therefore despite the large binding energy, 
this point is NOT a good candidate to focus the neutron 
scattering experiments. We will give a qualitative picture 
of wave-function in the following. 

Note the emergence of an instability at T— point near 
U c w 2.2t. Indeed at the level of a Hartree-Fock mean 
field one finds a transition to FM phase at 



1 



X °(q = 0,U = 0) 



2.231* 



However, real graphite for which U ~ it— 5* is stable. So 
the value of U c — 2.23* should be an artifact of RPA ap- 
proximation, and one expects by going beyond the RPA, 
to push U c above the 2.23* or equiyalently to obtain an 
screened value of U below 2.23* [lg]. Assuming that the 
renormalized value of U is less than 2.23*, we will do the 
rest of calculations for U — 2* ~ 5 eV. 



Once we find the location of resonance frequencies 
o; s (q), the next step is to calculate Z(q). The trick is 
to use the formula l|27|) . but with [^Rex°(q, ( q j 
given by direct differentiation of equation $5yp . that is 



8tt j 



[£k+ q + £k - ft^ s (q)] 2 



d 2 k 



V3U 2 a 2 \J [£ k +q + £k- ^(q)] 1 



(60) 



(61) 



Then using equations (|28|l and (|30|l the weight of neu- 
tron scattering peaks J(q), can be calculated, which is 
depicted in ngure|Hl The maximum value of the /(q) for 
q along TM, lies between ~ 0.1 rr = 0.2 YM w 0.3A" 1 
and - 0.2 rr = 0.4 YM « 0.61" 1 . At these wave vec- 
tors the energies of neutrons lie between ~ 0.5 eV and 
~ 1 eV which are much easier for neutron scattering than 
~ 2 eV at the M— point. The typical shape of the peak 
intensity along YK is similar to figure El Therefore ac- 
cording to this calculation, the best chance of detecting 
spin-1 collective mode, is at the points between YM and 
YK which are closer to Y, than K or M . The BZ bound- 
aries including K and M should be avoided. 

Note that for these points, I(q) becomes of the order of 
unity and the binding energy are typically Hu>b ^ 1 eV. 
The larger I(q) is, the easier will be exciting the triplet 
exciton. Equation l|70|l of appendix A, explicitly shows 
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FIG. 8: The approximate bound state wave-function within 
the RPA approximation. The center of mass momentum q 
is at the M— point. Unit of length is the lattice parameter 
a ~ 2.46A. Contours corresponding to |-0q(r) | > 0.1 have 
been plotted in the base. Note the anisotropic pattern of 
contours. The normalization of the wave-function is such that 
f/)q(0) = 1. At this point the spatial extent of the bound state 
wave-function is not so large and it can be identified with 
triplet excitons. This state can not be excited so easily. Since 
it requires momentum transfers beyond the first BZ, at which 
atomic form factors wash out the neutron peak. 



the relationship between spatial extent of the excitonic 
wave-function and Z(q). The small binding energies cor- 
respond to loosely bound particle-hole pairs and hence 
large wave-functions, which is according to l|7U|) synony- 
mous to small peak intensities. Intuitively speaking, it 
becomes harder to excite larger objects by neutrons of 
wave-length A ~ A. 

Note that very close to T— point, where in figure is 
not clear, according to (|54|l 7(q) sa Z(q) ~ q 3 which cor- 
responds to very large wave- functions. However as one 
expects form the qualitative behavior of binding ener- 
gies (figure and also the qualitative behavior of 7(q) 
(figure|H|), as q moves a bit away form T— point, the wave- 
function becomes of the order of a few lattice constants 
and hence visible by neutrons (figure 0). 

In figures[7|and|Hlthe contours show the region in which 
|^ q (r)| > 0.1 which gives a good feeling about the spa- 
tial extent of the bound-state wave-function. They re- 
spectively correspond to q = 0.2 x YM and q = TM. 
The normalization of the wave-function is such that 
ipq( r = 0) = 1- As can be seen form comparison of 
two figures, the spatial extent of the bound-state wave 
function for center of mass momentum corresponding to 
M— point, is of the order of two unit cells with highly 
anisotropic nature. 

Now let us discuss the effect of a bands. The a bands 



FIG. 9: \mx RP (w) as a function of lo at Q = Q.2TM » 
0.3A -1 . Note that u) in this figure is in units of t. The particle- 
hole continuum is clearly seen in this figure (dashed line). 
Below the continuum at w s (Q) » 0.30 t « 0.7 eV a very 
sharp peak in lmx RPA shows up (solid line). 



have their minima centered around the T— point. If we 
fit a quadratic dispersion to LDA data, we find that the 
inclusion of excitation form valence 7r band to a band do 
not modify the collective mode qualitatively. The height 
of the window below the particle-hole continuum of figure 
0] along the TK is not hight enough, so that the excita- 
tions to a do not shrink the window at all. However 
around M— point, at which the height of window is 2t, 
the inclusion of excitations to a band, reduce the height 
by maximum amount of ~ 0.5t at the M— point. But 
as can be seen from figure \5\ this does not open a decay 
channel for the normalized U ~ 2t. 



CONCLUSION 

We have evaluated the spin susceptibility of a graphcnc 
layer in RPA approximation for a short range interaction. 
This way we obtained a magnetic (spin-1) collective mode 
branch in non magnetic phase of graphite which exists in 
entire BZ with a very wide energy range form zero at T 
and if to ~ 2 eV near the M— point. This branch is 
below the particle-hole continuum and is protected from 
Landau damping. Therefore it might provide a mech- 
anism for transport of spin- only currents, over a wide 
energy range. We have shown that the long range tail 
of the Coulomb interaction of real graphite (3D stacks of 
graphene) does not destroy the conclusions drawn from 
the Hubbard interaction. We also presented the calcula- 
tion of the weight of the neutron peaks. The most ap- 
propriate region of the BZ to focus the neutron scatter- 
ing experiments is region in between T and M and region 
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mid-way between T and K— points (closer to T than M or 
K). Hot or epithermal neutrons (0.1 — 1.0 eV) are more 
appropriate for this region of momentum transfers. The 
dynamical form factor for a typical momentum transfer 
in the above region is schematically plotted in figure 
At the above mentioned regions, in one hand one does 
not require very high energy neutrons. On the other 
hand, the binding is strong enough which leads to small 
enough wave- functions which can be excited with hot or 
epithermal neutrons. 
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APPENDIX A: EXCITONIC WAVE-FUNCTION 

One can see more closely the relation between the Z(q) 
and the binding energies and real space profile of the 
wave-function as follows. At the RPA approximation, the 
eigenvector of the exciton state with spin S can be writ- 
ten as a linear combination of all the product eigenvectors 
of free charge carriers created in the crystal [12. fl3T | 



|q,S) £ 



k' 



<? q (k') |k'+q,-k';S) 



(62) 



where we have used the free electron- hole basis |k + 
q, — k; S) to expand the exciton wave- function. The 
triplet one which is relevant to our case is given by 



1 



|k + q,-k;l) = ^ (c^+qT^-kT + c k+qi d -ki 
+ ^ c k+qT«-kl + c k+q| a -kT 



10) 



where d_ k _ = c k(T creates a hole with spin projection 
—a and momentum — k. The exciton eigenvalue problem 
in free p-h basis becomes 

[£ c (k + q)-£"(k)-I¥ q ] ffq (k) 
+ 5>(k,k') 3q (k')=0 (63) 



Here E c (k + q) and E v (k) are total energies of single 
particle in conduction and valence band including their 
interaction energy with the other electrons of the band. 
The interaction between different pairs is included in in 
kernel C/(k, k') [l^. It can be seen from equation 164|) 
that in triplet channel kernel becomes a direct Coulomb 
term and moreover it is attractive. 

The RPA approximation amounts to ignore the self- 
energy effects and write (i.e. solving a two-body problem) 

£ c (k + q) -> e k+q/2 = e k+q/2 (65) 
£»(k)^4_ q/2 = -e k _ q/2 (66) 

and U(k, k') = —U for the Hubbard model. Note that 
wave-vectors are such that the wave-equation becomes 
manifestly time reversal invariant (k — > — k), as can be 
seen from the wave-equation below 

[Ek+q/2 + £k- q /2 - ftw s (q)] ffq(k) - U ^# q (k') = 

k' 

Let <7 q (k') = C, solve for <7 q (k) and sum over k to 
obtain the self-consistency equation 



C = 



uc 



£k+q/2 + ^k-q/2 ~ ^U>s (q) 



(67) 



which is exactly The RPA approximation, equation 143|l . 
The wave-function of the single particle-hole pair can be 
approximated by Fourier transform |isj of (? q (k). Thus 
we have 



^q(r) 



exp(ik.r) 



1 

N ^ £k+q/2 + £k-q/2 - ^s(q) 

A f rflv exp(ik.r) 

4tt 2 



£k+q/2 + &k-q/2 - 7^8 (<3L 



(68) 
(69) 



with A, being the unit cell area, which shows that the 
wave function is purely real, as it should be. Comparing 
equations (|60|l and (|68() shows that 

k J 

Note that we leave C in equation l|67|l is left undeter- 
mined. The normalization is such that ?Aq( r = 0) = 
1. Equation (|70|l shows that the intensity of collective 
mode peaks is a measure of inverse spatial extent of the 
particle-hole bound state. Smaller particle-hole bound 
states corresponds to sharper neutron scattering peaks. 



with kernel U(k, k') given by 

U{k,k') = -(ck + q;wk'|L/|ck' + q;wk) 

+2 ( 5 s ,o(ck + q;vk'|L/|vk;ck' + q) (64) 
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